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ABSTRACT 

In this paper we present the first results from an on-going serendipitous survey aiming 
to identify X-ray selected 'normal' galaxies (i.e. not AGN dominated) by combining 
archival XMM-Newton data with the Sloan Digital Sky Survey. In the first 4.5 deg 2 of 
this program we have identified a total of 11 'normal' galaxy candidates (8 of them 
with optical spectroscopy) with fluxes /jf(0.5 — 8keV) w 10~ 15 — 1CP 13 ergs -1 cm~ 2 . 
These sources are selected to have low X-ray-to-optical flux ratio (log fx / f op t ^ —2), 
soft X-ray spectral properties and optical spectra, when available, consistent with 
the presence of a stellar ionising continuum. These sources comprise both early and 
late type systems at redshifts z 0.2 with luminosities Lx(0.5 — 8keV) ~ 10 39 — 
10 42 ergs _1 . This dataset provides the first tight constraint on the surface density 
of X-ray selected 'normal' galaxies at relatively bright fluxes spanning two orders of 
magnitude (10~ 15 — 1CU 13 erg s -1 cm -2 ). The slope of the 'normal' galaxy log N — log S 
in the above flux range is estimated — 1.4±0.3 consistent with the euclidean prediction. 
We also discuss the prospects of 'normal' galaxy studies at X-ray wavelengths using 
both our continuously expanding survey and future X-ray missions. 
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1 INTRODUCTION 

In the local Universe 100 Mpc) 'normal' galaxies with 
X-ray emission dominated by stellar processes have been 
extensively studied using data from targeted observations 
of selected systems (e.g. Read, Ponman & Strickland 1997; 
Shapley, Fabbiano & Eskridge 2001; Zezas et al. 2003; Fab- 
biano et al. 2003; Jeltema et al. 2003). Although these in- 
dividual observations provide valuable information on the 
nature and the X-ray properties of 'normal' galaxies they 
cannot be used for statistical studies (e.g. log N — log S, lu- 
minosity function, clustering properties). 

It is only recently that ultra-deep Chandra surveys 
reaching fluxes /(0.5 — 2 keV) ~ 3 x 10~ 17 erg s -1 cm -2 (e.g. 
Brandt et al. 2001; Giacconi et al. 2002; Alexander et al. 
2003) have provided the first X-ray selected sample of 'nor- 
mal' galaxy candidates (Hornschemeier et al. 2003; Norman 
et al. 2004). These objects although too faint to signifi- 
cantly contribute to the X-ray background (XRB) are de- 
tected in increasing numbers with decreasing flux and are 
likely to outnumber powerful AGNs below /(0.5 — 2keV) ~ 
10~ 17 erg s -1 cm -2 (Hornschemeier et al. 2003). 

The studies above however, concentrate on 'normal' 
galaxies at faint X-ray limits (~ 10~ 16 erg s _1 cm -2 ) and 



moderate redshifts z = 0.1 — 1 (median ~ 0.3). At brighter 
fluxes (~ 10 -15 — 10 -13 erg s _1 cm -2 ) probing, on average, 
lower redshifts, z ^ 0.2, our knowledge on these systems 
remains scanty. This is because 'normal' galaxies at bright 
X-ray fluxes have low sky density. Therefore, wide area sur- 
veys are required to compile large enough samples to explore 
their statistical properties (luminosity function, surface den- 
sity; Georgakakis et al. 2004). Such a low- z study is essential 
to interpret the X-ray 'normal' galaxies in deeper X-ray sur- 
veys and to investigate their evolution. This is particularly 
important since the X-ray emission is the only observable 
window available that allows one to probe the evolution of 
the X-ray binaries and the hot gas in galaxies. This provides 
information on the star-formation history of the Universe 
that is complementary to that obtained from other wave- 
lengths (Ghosh & White 2001; Hornschemeier et al. 2002; 
Georgakakis et al. 2003; Norman et al. 2004). 

Wide areal coverage is difficult to achieve with the 
Chandra observatory due its limited field-of-view. On the 
contrary, XMM-Newton with 4 times larger field-of-view 
provides an ideal platform for such a study. In this paper 
we exploit the capabilities and the large volume of archival 
data of the XMM-Newton to serendipitously identify 'nor- 
mal' galaxy candidates in public fields selected to overlap 
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with the Sloan Digital Sky Survey (SDSS; York et al. 2000) 
to exploit the superb and uniform 5-band optical photom- 
etry and spectroscopy available in this area. In this paper 
we report the first results of this survey (Needles in the 
Haystack Survey; NHS) using a pilot sample covering a to- 
tal area of « 4.5 deg 2 in the flux range fx (0.5 — 8keV) ^ 
10~ 15 — 10~ 13 ergs -1 cm~ 2 . This dataset is continuously ex- 
panded as more SDSS fields and XMM-TVewfon archival 
data become available. Throughout this paper we adopt 
H = eSkms" 1 Mpc~\ Q M = 0.3 and fl A = 0.7. 



2 THE SAMPLE 

In this paper we use XMM-TVewfon archival observations, 
with a proprietary period that expired before September 
2003, that overlap with first data release of the SDSS (DRI ; 
Stoughton et al. 2002). Only observations that use the EPIC 
(European Photon Imaging Camera; Striider et al. 2001; 
Turner et al. 2001) cameras as the prime instrument oper- 
ated in full frame mode were employed. For fields observed 
more than once with the XMM- Newton we use the deeper 
of the multiple observations. A total of 20 fields are used 
with Galactic Nh in the range 1.4 — 13.2 x 10 20 cm~ 2 (mean 
3.7 x 10 20 cm~ 2 ) and PN good time intervals between 2 and 
67 ks (mean 18.5 ks). These observations are complemented 
with 8 XMM-Newton pointings from the northern region of 
the wide area, shallow XMM-A/e«rfon/2dF survey that also 
overlaps with the SDSS (Georgakakis et al. 2004). 

The X-ray data reduction, source detection and flux 
estimation are carried out using methods described in Geor- 
gakakis et al. (2004). Source extraction is performed in the 
0.5-8 keV merged PN+MOS images using the ewavelet 
task of SAS with a detection threshold of 5a. The extracted 
sources for each field were visually inspected and spurious 
detections clearly associated with CCD gaps, hot pixels or 
lying close to the edge of the field of view were removed. 
Figure 1 shows the solid angle covered by our survey as a 
function of the 0.5-8 keV limiting flux. 

To optically identify the X-ray sources we follow the 
method of Downes et al. (1986) as described in Georgakakis 
ct al. (2004) to calculate the probability, P, a given candi- 
date is the true identification. Here we apply an upper limit 
in the search radius, r < 7arcsec and a cutoff in the proba- 
bility, P < 0.05, to limit the optical identifications to those 
candidates that are least likely to be chance coincidences. 

'Normal' galaxy candidates are selected to have (i) ex- 
tended optical light profile (i.e. resolved) to avoid contami- 
nation of the sample by Galactic stars, (ii) X-ray-to-optical 
flux ratio log fx / fopt < —2, two orders of magnitude lower 
than typical AGNs (for the definition of log fx / fopt see 
Georgakakis et al. 2004) . The sample of 'normal' galaxy can- 
didates is presented in Table 1. 

The archival X-ray data used here include targeted ob- 
servations of nearby normal galaxies with low X-ray-to- 
optical flux ratio. Such sources have been excluded from 
Table 1. Moreover, a number of 'normal' galaxy candidates 
although not the prime target of the XMM-Newton pointing 
lie at the same redshift as the prime target and are there- 
fore most likely directly associated with it (e.g. cluster or 
group members). These sources are marked in Table 1 and 
are excluded from any statistical analysis. 
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Figure 1. Solid angle as a function of limiting flux desnity (5<r) 
in the total 0.5-8 keV band for our survey. 



3 OPTICAL AND X-RAY PROPERTIES 

Low luminosity Seyferts also lie in the log fx / fopt SS — 1 
regime and may contaminate our 'normal' galaxy sample 
(Lehmann et al. 2001). We explore this possibility using 
both the X-ray and the optical spectroscopic information 
available for our sources. 

For sources #2, 4, 6, 12 and 15 in Table 1 we can 
constrain the nature of the central engine using the diag- 
nostic diagrams proposed by Veilleux & Osterbrock (1987) 
that employ intensity line ratios that are least sensitive to 
dust reddening (H/3/[OIII] 5007, [S II] 6716+31/Ha, [OI] 
6300/Hq and [Nil] 6583/Ha). We adopt the classification 
scheme of Ho, Fillipenko & Sargent (1997) and find that 
sources #2, 6, 12 and 15 have HII type spectra while #4 is 
classified as a Seyfert 2. As will be discussed below the clas- 
sification of source #4 is consistent with its absorbed X-ray 
spectrum. 

For the remaining sources in Table 1 the classification 
using diagnostic diagrams cannot be applied because either 
they have absorption lines only (#9, 10, 16, 17), no spectro- 
scopic information (#1, 5, 7, 11, 14), or both the H/3 and 
the [OIII] 5007 A lines, essential for spectral classification, 
are absent from the optical spectrum (#3, 8, 13). The lat- 
ter sources have strong absorption features from a dominant 
evolved stellar population suggesting early and intermediate 
Hubble types. Galaxies with both absorption and emission 
lines have been extensively studied by Dressier et al. (1999) 
and Poggianti et al. (2000) on the basis of the equivalent 
widths of the H<5 (absorption) and [Oil] 3727A ( emission) 
lines. Our sources have EW~ —1 (absorption) and « 10 A 
(emission) for the H5 and [Oil] 3727A lines respectively. Fol- 
lowing the classification of Dressier et al. (1999) the evidence 
above suggests e(c)-type galaxies which are believed to be 
systems with constant star-formation rate over a period of 
time resulting in a mix of young and evolved stellar popu- 
lations. For these sources however, we cannot exclude the 
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Figure 2. Cumulative 'normal' galaxy counts. Filled circles are 
the 'normal' galaxy candidates in the present study. We plot 10 
instead of 11 points because sources #10 and 15 have the same 
flux and arc plotted as one point. The open circles are the source 
counts of Hornschemeier et al. (2003). The continuous line is the 
best fit to the log N — log S at bright fluxes from the present 
sample. The dashed line is the model prediction of the Georgan- 
topoulos ct al. (1999) X-ray luminosity function of HII galaxies 
assuming evolution of the form Lx ~ (1 + z) 3 . 



possibility of low-luminosity AGN activity diluted by the 
evolved stellar population (e.g. Severgnini et al. 2003) 

We can further constrain the nature of our 'normal' 
galaxy candidates via X-ray spectral analysis using the 
XSPEC vll.2 package. We use the C-statistic technique (Cash 
1979) specifically developed to extract information from 
even low signal-to-noise ratio spectra. The data are grouped 
to have at least one count per bin. An absorbed power-law 
model is fit to the data fixing the power law index to T = 1.9 
and leaving Nh free parameter. The results are shown in Ta- 
ble 1. The majority of the 'normal' galaxy candidates have 
soft X-ray spectral properties suggesting obscuration consis- 
tent with the Galactic Nh- The only exception is source #4 
which is significantly absorbed. 

The final sample of 'normal' galaxy candidates that 
are not associated with the prime target of a given XMM- 
Newton pointing and not showing evidence (e.g. X-ray or 
optical spectra) for obscured AGN activity comprises 11 
sources, #1, 2, 3, 5, 6, 8, 9, 10, 11, 13, 15 in Table 1. Of 
these 11 sources only sources #1,5 and 11 do not have op- 
tical spectroscopy available. 



4 DISCUSSION 

In this paper we report on an on-going program aiming to 
identify X-ray selected 'normal' galaxies in the flux range 
/(0.5 - 8keV) w 10~ 15 - 10~ 13 ergs" 1 cm" 2 by combin- 
ing archival XMM-Newton data with the SDSS. At present 
our on-going survey covers an area of 4.5 deg 2 and has 
identified a total of 11 'normal' galaxy candidates with 
log fx/fopt < —2. Unlike deeper samples our sources are rel- 
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Figure 3. Lx(0.5 — 8kcV) against redshift. Filled circles are for 
our sample of 'normal' galaxy candidates (including sources asso- 
ciated with the prime target of a given XMM-Newton pointing). 
Open circles represent the sample of Hornschemeier ct al. (2003). 
The two surveys cover different but complementary regions of the 
Lx — z space. 



atively bright at both the optical and the X-ray wavelengths 
allowing detailed study of their nature thus, minimising any 
AGN contamination. Both the soft X-ray spectra and the 
optical spectroscopic properties are consistent with the ab- 
sence of nuclear activity. 

Figure 2 plots the log N — log S of our sample in com- 
parison with that of 'normal' galaxy candidates identified 
in the Chandra Deep Field North survey probing fluxes w 



10 



~ 16 erg s 1 cm 



(Hornschemeier et al. 2003). The present 
study provides the only constraints on the number density 
of 'normal' galaxies at relatively bright fluxes in the range 
/(0.5 - 8keV) « 10" 15 - 10~ 13 ergs" 1 cm" 2 spanning two 
orders of magnitude. The differential 'normal' galaxy counts 
in the 0.5-8 keV band are fit by a power law yielding a slope 
of — 1.4 ± 0.3 for the cumulative number counts. This is 
in excellent agreement with the results of Hornschemeier 
et al. (2003) at fainter fluxes (-1M±%H). Figure 2 also 
shows that our wide area medium-deep serendipitous sur- 
vey is complementary to the ultra-deep pencil beam Chan- 
dra observations. This is further demonstrated in Figure 3 
plotting the L — z plane for both the present sample and 
that of Hornschemeier et al. (2003). The two surveys cover 
different but complementary regions of the parameter space 
with our sample probing the low-2, high-Lx regime. Such 
a wide L — z plane coverage is essential to explore the evolu- 
tion of 'normal' galaxies at X-ray wavelengths. Although our 
sample is not large enough yet for such a study, we plan to 
expand it to ~ 20 deg 2 by exploiting the 2nd data release of 
the SDSS. Our enlarged sample will comprise fa 40 'normal' 
galaxies sufficient for luminosity function estimates provid- 
ing an anchor point at low redshifts, essential for comparison 
with higher-2 studies (e.g. Hornschemeier et al. 2003; Geor- 
gakakis et al. 2003; Norman et al. 2004). 

However, even the deepest X-ray samples todate find 



© 0000 RAS, MNRAS 000, 000-000 



4 Georgakakis et al. 



ID 


a x 


Sx 


P 


r 


z 


net 




(J2000) 


(J2000) 


(per cent) 


(mag) 




counts 


l b 


03 58 05.22 


+01 09 48.96 


0.01 


14.93 


0.074 


53 


2 


08 30 59.21 


+52 37 47.75 


0.80 


17.69 


0.1363 


48 


3 


08 31 14.54 


+52 42 24.77 


0.01 


15.24 


0.0641 


120 


4 


08 31 39.09 


+52 42 05.51 


0.02 


15.70 


0.0587 


337 



/x(0.5-8kcV) 

- 14 n r rr o 1 2 \ 



L x (0.5-8keV) N H type a 

(10 41 ergs" 1 ) (10 22 cm~ 2 ) 



1.52 + 0.34 


2.36 


< 0.04 


early 


0.34 + 0.13 


2.01 


< 0.13 


NL 


0.56 + 0.10 


0.66 


< 0.02 


EA 


2.07 + 0.14 


2.00 


14 9+ 09 


NL 



5 b 


08 


32 


02.58 


+52 


47 


13.74 


0.01 


16.01 


0.105 


55 


6 


08 


32 


28.14 


+52 


36 


22.36 


<0.01 


13.82 


0.0170 


716 


>jb,c 


09 


17 


58.59 


+51 


51 


04.75 


0.31 


16.66 


0.219 


22 


8 


09 


35 


18.39 


+61 


28 


31.58 


0.27 


16.50 


0.1241 


54 



0.25 + 0.09 


0.81 


< 0.1 


late 


6.43 + 0.29 


0.52 


< 0.02 


NL 


0.61 + 0.26 


9.76 


< 0.06 


early 


0.96 + 0.25 


4.48 


< 0.08 


EA 



9 


09 36 19.46 


+61 27 23.11 


0.06 


16.55 


0.1311 


67 


0.95 + 0.21 


5.25 


< 0.03 


AB 


10 


12 31 47.21 


+64 14 01.68 


< 0.01 


13.01 


0.009 


67 


1.95 + 0.14 


0.047 


< 0.02 


AB 


11» 


12 32 53.15 


+64 08 50.93 


0.17 


15.43 


0.139 


50 


1.03 + 0.26 


6.12 


< 0.06 


late 


12 c 


12 45 31.99 


-00 32 08.59 


0.01 


13.12 


0.0055 


353 


2.83 + 0.19 


0.02 


11+ - 07 
u - li -0.03 


NL 



13 


12 44 52.26 


-00 25 49.66 


0.01 


15.34 


0.0823 


95 


0.57 + 0.12 


1.10 


< 0.03 


EA 


14&,c 


13 03 01.13 


+67 25 17.76 


0.16 


16.16 


0.109 


49 


1.38 + 0.35 


4.90 


001 -o.oi 


early 


15 


13 41 37.85 


-00 25 55.30 


0.05 


16.86 


0.053 


30 


1.95 + 0.52 


1.44 


< 0.09 


NL 


16 c 


23 53 40.52 


-10 24 19.91 


0.02 


15.07 


0.0741 


52 


2.52 + 0.52 


3.92 


< 0.03 


AB 


17 c 


23 54 05.69 


-10 18 30.64 


0.06 


15.70 


0.0734 


50 


1.06 + 0.33 


1.61 


< 0.02 


AB 


"AB: 


absorption lines; NL: Narrow 


emission lines: EA: both 


narrow 


emission 


and absorption lines 









^photometric redshift and classification from the SDSS 

c source at the same redshift as the target of the XMM- Newton pointing. 

The columns are: 1: identification number; 2: right ascension of the X-ray source; 3: declination of the X-ray source; 

4: probability the optical counterpart is a spurious alignment; 5: optical magnitude; 6: spectroscopic or photometric redshift from the SDSS; 
7: 0.5-8 keV net counts; 8: X-ray flux in the 0.5-8 keV spectral band in units of 10 — 14 ergs -1 cm~ 2 ; 

9: 0.5-8 keV X-ray luminosity in units of 10 41 ergs _1 ;10: Njj from the X-ray spectrum assuming V = 1.9; 11: spectral type. 



Table 1. The candidate 'normal' galaxy sample. 



'normal' galaxies at relatively low-z (median redshift z ~ 
0.3; Hornschemeier et al. 2003) significantly lower than the 
star-formation peak at z « 1 — 2. This underlines the need 
for deeper X-ray observations to use X-ray selected 'normal' 
galaxies as cosmological probes. The ESA XEUS mission 
with its unparalleled sensitivity will provide a great leap for- 
ward in these studies. Assuming the goal resolution of about 
2 arcsec FWHM for its mirrors, XEUS will attain, in just 
20 ks, about equal sensitivity to the deepest (2 Ms) Chandra 
exposure. Extrapolation of the log N — log S relation derived 
here to fx(0.5 — 8) = 10~ 16 ergcm~ 2 s _1 (i.e. approximate 
flux limit of the 2 Ms Chandra Deep Field North; Alexander 
et al. 2003) , suggests that about 10 galaxies will be detected 
in the 25 arcmin 2 XEUS Wide Field Imager. Therefore, a few 
tens of typical exposures with XEUS will readily accumulate 
a large sample of X-ray selected galaxies for the study of the 
galaxy luminosity function and evolution at redshifts z ^ 1 
with a median redshift of z ~ 0.3. Deeper exposures are nec- 
essary to probe the star-formation peak at higher redshifts. 
In a deep 200 ks exposure, XEUS is expected to probe fluxes 
as faint as /x(0.5 - 8keV) « 5 x 10~ 18 ergs -1 cm~ 2 and to 
detect more than 500 galaxies (an estimate obtained by ex- 
trapolating our log N — log S). At these flux levels the mean 
redshifts probed are z ~ 1 as derived from the X-ray lumi- 
nosity function of Georgantopoulos et al. (1999) assuming 
evolution of the form Lx ~ (l + z) a . Therefore, combination 
of the local XMM- Newton galaxy samples with the XEUS 
data will allow an unprecedented study of the star-formation 
history at X-ray wavelengths. 
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